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Abstract 

We present a new parametrization of the minimal seesaw model, expressing the heavy- 
singlet neutrino Dirac Yukawa couplings {Y v )^ and Majorana masses M^. in terms of ef- 
fective light-neutrino observables and an auxiliary Hermitian matrix H. In the minimal 
supersymmetric version of the seesaw model, the latter can be related directly to other low- 
energy observables, including processes that violate charged lepton flavour and CP. This 
parametrization enables one to respect the stringent constraints on muon-number violation 
while studying the possible ranges for other observables by scanning over the allowed pa- 
rameter space of the model. Conversely, if any of the lepton-flavour-violating process is 
observed, this measurement can be used directly to constrain (Y v )ij and M Ni . As applica- 
tions, we study flavour-violating r decays and the electric dipole moments of leptons in the 
minimal supersymmetric seesaw model. 
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1 Introduction 



Experiments on both atmospheric and solar neutrinos have now provided 'smoking guns' for 
neutrino oscillations. Most recently, the direct SNO measurement of the solar neutrino flux 
via neutral-current scattering confirms solar neutrino oscillations Q and favours strongly 
the LMA solution 0. This region of parameter space is within reach of the KamLAND 
experiment, and is expected to be probed soon ||. The existence of large mixing angles for 
both solar and atmospheric neutrinos J| is one of the biggest mysteries in particle physics. 
The most favoured mechanism for generating neutrino masses is the seesaw mechanism ||, 
which naturally explains their small sizes. However, it is an open question whether the seesaw 
mechanism can explain why mixing in the lepton sector seems to be larger than in the quark 
sector. In the absence of a theory of flavour, it is important to study the consequences of 
neutrino mixing for as many physical observables as possible. 

In the minimal supersymmetric seesaw model, lepton-flavor-violating (LFV) phenomena 
provide a tool to study indirectly neutrino parameters and probe other aspects beyond the 
large mixing angles measured in neutrino oscillations. If supersymmetry breaking originates 
from physics beyond the heavy singlet neutrino mass scale, LFV slepton masses are in- 
duced radiatively |], [7| via the Dirac Yukawa couplings of the neutrinos, even if the input 
supersymmetry-breaking parameters are flavor-blind. On the other hand, the light neutrino 
masses and mixings depend on both the Yukawa couplings and the Majorana masses of the 
heavy singlet neutrinos. Thus one can hope to reconstruct the physical parameters in the 
heavy singlet-neutrino sector entirely in terms of the light neutrino data and low-energy 
observables such as rates for LFV processes |§. To this end, in this paper we present a 
parameterization of the minimal seesaw model and apply it to the minimal supersymmetric 
version of the seesaw model. 

The essence of our parametrization is the following. The minimal seesaw mechanism, 
whether supersymmetric or non-supersymmetric, involves 18 physical degrees of freedom, 
including 6 real mixing angles and 6 CP-violating phases. On the other hand, the induced 
light-neutrino mass matrix has 9 degrees of freedom, including 3 real mixing angles and 3 CP- 
violating phases. Thus we need 9 additional degrees of freedom to parametrize completely the 
seesaw mechanism. These can be chosen in such a way as to be related to low-energy leptonic 
observables in the supersymmetric version of the seesaw model. We recall that the LFV 
renormalization of the supersymmetry-breaking parameters at low energy are proportional 
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to 

Hij = ETO«(n)«log-^-, (1) 

k M N k 

where (Y u ) and Mn are the heavy singlet-neutrino Dirac Yukawa couplings and Majo- 
rana masses, respectively, and Mq is the GUT scale where the initial conditions for the 
supersymmetry-breaking parameters are imposed. Since if is a Hermitian matrix, it has 
9 degrees of freedom including 3 real mixings and 3 phases. This implies that we can 
parametrize the seesaw mechanism by the light neutrino mass matrix M. v and the Hermi- 
tian matrix H according to 

(M V ,H)—>(Y V ,M N ). (2) 

As a result, we can obtain Y v and Mn that yield automatically the light neutrino masses 
and mixings measured in oscillation experiments. However, the main motivation for our 
parametrization comes from its power in studies of the charged-lepton physics in the super- 
symmetric seesaw model []. 

The LMA solution to the solar neutrino anomaly tends to predict a large branching 
ratios for fi — > ej in the supersymmetric seesaw model |9|, [□]], which may be within 
reach of near-future experiments, or even beyond the current experimental bound Q. This 
does not imply that the supersymmetric seesaw model is strongly constrained, because it 
has a multi-dimensional parameter space. However, it is difficult to scan efficiently over 
the allowed parameter space while satisfying the /i — > ey constraint. Our parametrization 
solves this difficulty, because the parameter matrix H is related to the solutions of the 
renormalizat ion-group equations. Therefore, it is straightforward to choose a parameter 
region where /i — > is suppressed, but the other low-energy observables may vary over their 
full ranges. Furthermore, if some future experiment discovers a LFV process or the electric 
dipole moment (EDM) of some lepton, this observation will be directly related to H and 
thus to the neutrino parameters. In our parametrization, the high-energy neutrino couplings 
and masses can be expressed entirely in terms of the induced low-energy observables. 

Our work is organized as follows. In Section 2 we outline the new parametrization and 
our procedure for analyzing charged-lepton decays. In Section 3 we explain the relation 
between our parametrization and the physical observables. In Section 4 we present a study 

*We emphasize, though, that the parametrization itself is more general, and does not depend on the 
existence of supersymmetry. 

2 Also, some explicit models predict the third neutrino mixing parameter U e 3 to be O(10 _ ^ 1_2 - ) ), which 
may also lead to a large branching ratio for /i — > ej |^2| . 
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of LFV r decays and the EDMs of the electron and muon in the supersymmetric seesaw 
model, as applications of our approach. We find that r — > fJ,(e)"f can saturate the current 
experimental bound, even when /i — > is suppressed enough to be acceptable. The EDMs 
of the muon and electron generally fall below 10~ 27 (10 -29 )e cm in our random parameter 
scan. We also present the relation between Br(r — > /i(e)7) and Br(r — > fi(e)f + fr). Section 
5 summarizes our conclusions. 



2 Parametrization of Neutrino Couplings and Masses 

In view of the subsequent application to the supersymmetric version of the seesaw model, 
we illustrate the parametrization for this case, though it is also valid in the absence of 
super symmetry. The leptonic superpotential of the supersymmetric version of the minimal 
seesaw model is 

W = N^(Y l/ ) ij LjH 2 + E^(Y e ) ij LjH 1 + ^(M^N? , (3) 

where the indices i,j run over three generations and (M^)^ is the heavy singlet-neutrino 
mass matrix. In addition to the three charged-lepton masses, this superpotential has 18 
physical parameters, including 6 real mixing angles and 6 CP-violating phases. 

At low energies, the effective superpotential obtained by integrating out the heavy neu- 
trinos is 

W eS = E^Y^LjH, + 1 (MMLiH 2 )(LjH 2 ) , (4) 

2v z sin p 

where we work in a basis in which the charged-lepton Yukawa couplings are diagonal. The 
second term in (^) leads to the light neutrino masses and mixings. The explicit form of A4 U 
is given by 

( KA \ \ - (Yv(Qk))ki{Yy{Qk))kj 

k M N k 

where the heavy-singlet neutrino Dirac Yukawa couplings Y v and masses Mat. are defined 
at the renormalization scale Qk = M^ kl and in our notation < Mn 2 < Mn 3 . It is 
important to distinguish between the renormalization scales for different components in the 
Yukawa coupling matrix, since the EDMs of charged leptons in the supersymmetric seesaw 
model are sensitive to non-universal radiative corrections to the supersymmetry-breaking 
parameters, which come from the non-degeneracy of the heavy singlet neutrino masses ||13|| . 
For simplicity, we ignore the renormalization of M. u after the decoupling of the singlet 
neutrinos. 
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The light neutrino mass matrix M. v (|5|) is symmetric, with 9 parameters, including 3 real 
mixing angles and 3 CP-violating phases. It can be diagonalized by a unitary matrix U as 

U T M V U = M° . (6) 

By redefinition of fields one can rewrite U = VP, where P = diag(e J<?il , e*^ 2 , 1) and V is the 
MNS matrix, with the 3 real mixing angles and the remaining CP- violating phase. 



The key proposal of this paper is to characterize the seesaw neutrino sector by M. v and 
a Hermitian matrix H, whose diagonal terms are real and positive, which is defined in terms 
of Y v and the heavy neutrino masses by 

Ha = Y,(V{Qk))ki{Y v {Q k )) kj \og^-, (7) 

k M N k 

with Mq the GUT scale. The Hermitian matrix H has 9 parameters including 3 phases, 
which are clearly independent of the parameters in M. v . Thus M. v and H together provide 
the required 18 parameters, including 6 CP-violating phases. 

Although our parametrization also includes an unphysical region, it has the merit of 
suitability for comprehensive studies of the minimal supersymmetric seesaw model. In this 
model, the non-universal elements in the left-handed slepton mass matrix, which induce 
the charged LFV observables, are approximately proportional to H if the slepton masses 
are flavour independent at Mq. Thus, this parameterization allows us to control the LFV 
processes and scan over the allowed parameter space at the same time. Conversely, if some 
LFV process is discovered in the future, its measurement can be incoprorated directly into 
our parametrization of the neutrino sector. 

We now explain how to reconstruct the heavy singlet-neutrino sector from knowledge of 
A4v and H. First we recall the parametrization of the neutrino Dirac Yukawa coupling given 
in H, 



(Y*(Qi)) 



v sin P 
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where R is an auxiliary complex orthogonal matrix: RR T = R T R = 1. Using this parametriza- 
tion, H becomes 

1 



H = 2 • 2a UJM v tfM N RjM v U^ (9) 
v 2 svn P v v 



where Mjq. = M^. log(Mc/Mjv.). If we can diagonalize the following Hermitian matrix H', 



— i / 1 



H' = JM V WHUJMu v 2 sm 2 (3, (10) 
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by the complex orthogonal matrix R': 



H' = R M^R', 



(11) 



then we can calculate the heavy singlet neutrino masses from Mn and the corresponding Y v 
from @ taking R = R '. 

However, the Hermitian matrix H' cannot always be diagonalized by a complex orthog- 
onal matrix: the condition for such a diagonalization is that all the eigenvalues of H'*H' are 
positive, in which case R' is given by the eigenvectors of H'*H'. This reflects the fact that 
our parametrization also includes an unphysical region, so that every chosen H does not 
necessarily give physical neutrino masses and couplings. Since our objective in this paper is 
to survey the multi-dimensional parameter space using scatter plots, this shortcoming is not 
critical. 

In our subsequent analysis, we first generate randomly the matrix H, the phases and 
the common mass scale in the light neutrino sector, and then calculate the correspond- 
ing heavy neutrino masses and couplings. The Yukawa couplings (Y v )ij contribute to the 
renormalizat ion-group (RG) equations above Mjv j? since the corresponding singlet neutrino 
is dynamical there. When we derive the Yukawa couplings at the GUT scale, we introduce 
(Y„)ij in the RG equations at Qi = Mjy. where the neutrinos appear. When evaluating 
the supersymmetry-breaking parameters at the weak scale, the right-handed neutrinos are 
integrated out at their own mass scales. 

3 Observables 

In the previous Section we presented our parametrization of the minimal seesaw mechanism 
in terms of the light-neutrino mass matrix M. v and a Hermitian parameter matrix H. Here we 
make explicit the correspondence between this parametrization and low-energy observables 
in the super symmetric version of the seesaw model. 

3.1 Neutrino Experiments 

As already mentioned, the light-neutrino mass matrix Ai u contains nine physical parameters: 
3 mass eigenvalues, 3 mixing angles, 1 CP-violating mixing phase in the MNS matrix, and 
2 CP-violating Majorana phases, the LMA solution is now favoured, following the SNO 
neutral-current result. Thus, the favoured regions for the atmospheric and solar neutrino 
parameters are 



Am% 



32 




(12) 



5 



sin 2 20. 



23 



(0.8-1.0) 



Am 2 21 



tan 2 6» 



: l -( 4 - 5 )eV 2 
l2 - (0.2-0.6). 



The CHOOZ |14j and Palo Verde |15] experiments provide the constraint 



sin 2 2^13 < 0.1. 



(13) 
(14) 
(15) 



(16) 



These parameters, together with the CP-violating mixing phase in the MNS matrix, may 
be measured in future experiments, such as the KamLAND and the neutrino factory. There 
would still be three undetermined parameters, the normalization of the neutrino mass and 
the Majorana phases. The neutrinoless double beta decay matrix element is proportional to 



(17) 



and so would provide a constraint on the neutrino mass scale and Majorana phases, if it 
could be measured. 



3.2 Charged LFV Processes 

If the supersymmetry-breaking parameters at the GUT scale are universal, off-diagonal com- 
ponents in the left-handed slepton mass matrix and the trilinear slepton coupling A e are 
induced by renormalization, taking the approximate forms 



?V. 

L> % 3 



(5m 



1 



i7l z 
1 



(3m 2 + A 2 )K, 



y • 



57T 



-AoY ei Hij , 



(18) 



where i ^ j, and the off-diagonal components of the right-handed slepton mass matrix are 
suppressed. The parameters m and A are the universal scalar mass and trilinear coupling 
at the GUT scale. Here, we ignore terms of higher order in Y e , assuming that tan/? is not 
extremely large. Thus, the parameters in H may in principle be determined by the LFV 



processes of charged leptons. Currently, /i 
on them: 



ej experiments give the following constraints 



H 
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< 



10~ 2 x tan"2 (3 



10 _1 x tan"2 @ 



m \ 2 / Br(/i — > ej) 
lOOGeV/ V 1.2 x 10- 11 

m \ 2 / Br(/x — > ej) 
lOOGeV/ I 1.2 x 10" 11 



(19) 
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where we take (mj) 



m§. These components may also be measured directly in future 



collider experiments, if the sleptons are produced there |16], [17] . 

Although the matrix H has three CP-violating phases, two of them are almost irrelevant 
to charged LFV phenomena. The two phases may be moved from H to A4 V by a rotation of 
L. In fact, there is only a single Jarlskog invariant obtainable from H 



J 



lmHi 2 H23H 3 i 



(20) 



which determines the T-odd asymmetry in fi — ► 3e (JTJ 

We kept in ([Of) only the leading-order contributions to the soft supersymmetry-breaking 
parameters, and ignored higher-order corrections. If some components of H are suppressed, 
non-trivial flavour structure may emerge in the higher-order corrections. At O (log 2 Mg/M^ 3 ) 



or (9 ( log Mq I M jv 3 togM^/AfjvJ (i 7^ j) 

1 



(ml 



and (A e ) get the following corrections: 



(S'm 



(<5'A 



(4tt) 4 ' 

1 



(21) 



(4vr) 4 

k<l 



A Y ei 



6 {X fc , X,} log ^ + 4 [X fc , X] log ^ 



log 



M, 



G 



M: 



N 3 



where 



ij 



(Y;(M Nk )) ki (Y u (M Nk )) kj . 



(22) 



(23) 



The second term in ([H]) and the term in (|22]) come from threshold corrections at the heavy 
singlet-neutrino scale, due to the non- degeneracy of their neutrino masses. 



3.3 EDMs of the Charged Leptons 

The EDMs of the charged leptons depend non-trivially on the parameters in the supersym- 
metric seesaw model. If the heavy singlet-neutrino masses are degenerate, the EDMs of the 
charged leptons are strongly suppressed by the chiral structure of the seesaw model, and are 
proportional to 



Im 



yy^y fyty y^y]y^y 



(24) 



This is similar to the neutron EDM in the Standard Model. However, when the heavy 
singlet-neutrino masses are not degenerate, the EDMs may be enhanced significantly. The 
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trilinear coupling A e gets a threshold correction at the heavy singlet neutrino scale, and may 
get radiatively-induced diagonal phases proportional to 



ImiX^X^u log M N JM N . ^ 0, 



as m 



13| . This depends non-trivially on the CP- violating phases, including the two 



Majorana phases in M. v and two phases in # that are irrelevant for LFV. 



4 Phenomenological Analysis 

Using the parametrization proposed above, we now study the branching ratios for LFV r 
decays, such as Br(r — > /ry), Br(r — ■> ej) and Br(r — > 3/), and the EDMs and d e . We fix 
the light neutrino parameters: Am^ = 3 x 1CT 3 eV 2 , Am^ = 4.5 x 1CT 5 eV 2 , tan 2 #23 = 1 
and tan 2 # 12 = 0.4 corresponding to the LMA solution for the solar neutrino anomaly. Since 
the bound on the angle #13 is quite stringent, our results depend very weakly on its actual 
value. We fix sin #13 = 0.1 and 5 = n/2. We study both the normal and the inverse hierarchy 
of light neutrino masses, since neutrino oscillations do not discriminate between these two 
cases As input parameters, we then have the lightest effective neutrino mass mi (or m 3 
for inversely ordered neutrinos), which we generate in the range (10~ 4 — 0.3) eV, the two low 
scale Majorana phases 01,2 and the matrix #, which we generate randomly. 
We study two different limits of the parameter matrix #, of the form 




ifi= b d , (25) 



and 




#2 = 6 0, (26) 



where a, b, c are real and positive, and d is a complex number. We sample these parameters 
randomly in the range 10~ 2 < a, b, c, \d\ < 10, with distributions that are flat on a logarithmic 
scale. Also, we require the Yukawa coupling- squared to be smaller than Air, so that Y v remains 
perturbative up to M G . 

In the above ansatz, we take #12 = and H 13 H 32 = because these conditions suppress 
Br(/i — > ej), as seen in (0). We show Br(/x — > ej) as a function of the heaviest right-handed 
3 Future neutrinoless double-beta decay and oscillation experiments will resolve this ambiguity. 
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neutrino mass Mn 3 for the two structures Hi and H 2 in Fig. [I]. As a reference, here we take 
m 1/2 = 300 GeV, m = 100 GeV, A = -300 GeV, tan/5 = 10 and sign(fi) = +1. 

With the chosen forms H\ and H 2 , Br(/i — > cy) is suppressed in a broad range of param- 
eters. In these figures /i — > ey is induced entirely by C(log M G /M N3 log M N JM N .) {i ^ j) 
corrections to the slepton mass matrix and the slepton trilinear coupling . Thus, if all 

non-zero components in H\ or H 2 are of order unity, the correction to (m|)i2 is not necessary 
negligible, and one may find values of the branching ratio above the present experimental 
bound. We also see in Fig. that improving the present sensitivity to Br(/i — > ej) by three 
orders of magnitude, which is currently being undertaken at PSI [ 19| , would be interesting 
for a large fraction of the models studied. The MECO and PRISM searches for /i~ — e~ 
conversion on nuclei pO] would also be interesting in this respect, and the sensitivities to 



both processes could be improved at the front end of a neutrino factory [21 




Figure 1: Scatter plot of Br(fi —>■ ej) against the heaviest singlet neutrino mass for the 
ansatz (a) H x and (b) H 2 . We take m 1/2 = 300 GeV, m = 100 GeV, A = -300 GeV, 
tan/3 = 10 and sign(^) = +1. Other input parameters are specified in the text. 



The ansatz Hi minimizes r — > cy while r — > /i7 can be large, and the opposite is the 
case for H 2 , since these processes are sensitive to Hi 3 and H23, respectively. Similarly, the 
EDM of the muon can be maximized in Hi while the electron EDM can be large in H 2 . We 
exhibit these results in the following subsections. 
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4.1 LFV t Decay 

In Fig. we present Br(r — > ^7) for the ansatz Hi, assuming either the normal or the inverted 
hierarchy for the light neutrino mass spectrum. The horizontal axis is the lightest stau mass 
m fl , and the other supersymmetry-breaking parameters are determined by choosing the 
SU(2) gaugino mass to be 200 GeV, A Q = 0, \i > 0, and tan/? = 10 and 30. The parameters 
in M. v and H are the same as in Fig. p]. The branching ratio scales as tan 2 /?. We see from 
these figures that the branching ratio is similar for the normal and inverted hierarchies of 
light-neutrino masses. In our parametrization the branching ratio is determined mainly by 
H and the sparticle mass spectrum. The dependence on the details of M. v appears through 
the heavy singlet neutrino masses, which influence the branching ratios only logarithmically. 
We find that Br(r — > ^7) can reach even above the present experimental bound, attaining 
10" 4 (10~ 5 ) for tan/3 = 30(10), these limits are arising from the perturbative bound on the 
neutrino Yukawa coupling. 

Next we show Br(r — >• cy) for the ansatz H2 in Fig. ^. The input parameters are the 
same as in Fig. 0, and the behaviour is similar to that of Br(r — > fij) in Fig. 0. This process 
may also reach to 10~ 4 (10~ 5 ) for tan/3 = 30(10). 

We see in Fig. ^ § that improving the present sensitivity to Br(r — > (t/ej) by two 
orders of magnitude, which seems feasible with the present generation of B factories |p2| . 
would be very interesting for many of the models studied. Whilst the B factories have good 
sensitivities to r — > /I'y and r — > ej, hadron colliders may have better reach for r — > 3£ 
decays. We present here a simple formula for the LFV r decays to three charged leptons. 
When (m|) 2 3 is non-vanishing, non-zero r — > fi'j and r — > fi£ + £~ are both predicted. The 
photonic penguin diagram tends to dominate over other contributions in the trilepton final 
state due to the phase-space integral. When tan (3 is large, the dominance is even stronger. 
When the photonic penguin diagram is dominant in r — > fi£ + £~ , 

B < r - 3 ^ = fL5(iogl£_iiUJL (27) 
Br(r -> fij) 8tt 3 \ m 2 4 J 440 ' v ; 

B ^ - /^) = £? L^jLi (28) 
Br(r -> ^7) 8tt3\ & m 2 e 3 J ~ 94 ' 1 ' 

The branching ratio for fie + e~ is larger than that to fifi + [i~, because the phase space is 
larger. Similarly, we get the following relations for the r-e transition: 

Br(r -> 3e) 1 



Br(r -> erf) ~ 95 ' 
Br(r -> e2fi) _ 1 
Br(r ->• e7) ~ 430 



(29) 
(30) 



10 



10" 4 



10" 4 



tan|3=30 (normal hierarchy) 




m~ [GeV] m~ [GeV] 



Figure 2: Scatter plot of Br{r — > /ry) against the lightest stau mass for the ansatz Hi. We 
take the SU(2) gaugino mass to be 200 GeV, A = 0, ji > 0, and tan/3 = 10 and 30. We 
consider both the normal and inverted hierarchies for the light neutrino mass spectrum. 
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Figure 3: Scatter plot of Br{r — > ej) against the lightest stau mass for the ansatz H 2 . The 
input parameters for the supersymmetry-breaking parameters are the same as in Fig. f|. 
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Thus, from the experimental bound on Br(r — > fi/ej), we find that Br(r — > /x2e) and 
Br(r — > 3e) can reach 10~ 8 , and Br(r — ► e2/i) and Br(r — > can reach 10~ 9 . 

4.2 EDMs of Charged Leptons 

The EDMs of charged leptons are induced essentially through the threshold correction to A e 
at the heavy singlet-neutrino scale, and the dependence on H and A4 U is complicated. From 
(f22|) it is found that the imaginary parts of diagonal terms in A e can reach O(0. 1)% if the 
(Xk)ij are no larger than 4ir. In this case, the muon and electron EDMs can roughly reach 
the level of 10 _25 e cm and 10 _27 e cm, respectively. However, the non-negligible contribution 
of the heavy singlet neutrino threshold correction may lead to large Br(/i — > ej) even if H\2 
and Hi 3 H 3 2 are suppressed. Therefore careful numerical study is required for the predictions 
of EDMs in the super symmetric seesaw model. 

We show first in Fig. [| predictions for the EDMs of the muon and electron from a 
random sampling of the parameter space, as a scatter plot of the muon EDM against the 
left-smuon mass for the ansatz Hi. As before, we take the SU(2) gaugino mass to be 
200 GeV, A = — 3m , \x > 0, and tan/3 = 10. The parameters in M. v and H are the 
same as in Fig. [l| We assume the normal hierarchy for the light neutrino mass spectrum, 
and impose the constraints from Br(/x — > ej) and Br(r — > /x(e)7). We find that the muon 
EDM can reach lO - ^ 27-28 ^ cm in this sampling, which is reasonable, since (m|)i2 should 
suppressed by 10 - ( 3-4 ) because of the experimental bound on Br(/i — > ej). The EDM 
increases proportionally to A , as expected, but is insensitive to tan/3, if it is not extremely 
large. 

Similarly, in Fig. |5| we show the electron EDM plotted against the left-selectron mass for 
the ansatz H2. We find that the electron EDM can reach lO^ 29-30 -^ cm in this sampling. 
Again, the main limiting factor is the experimental bound on Br(/i — > c-f). 

We recall that a proposal has been made to BNL that aims at a sensitivity of 10~ 24 e cm 
for the muon EDM [^], and the front end of a neutrino factory may be able to reach a 
sensitivity of 5 x 10~ 26 e cm [U\. On the other hand, although the present upper limit on the 



electron EDM is 1.6 x 10 27 e cm [E4[], a technique has been proposed that may be sensitive 



to 10 32 e cm p5| . This would be sensitive to many of the models studied. 



5 Conclusions 

We have presented a new parametrization of the minimal seesaw mechanism, which enables 
the heavy neutrino Dirac Yukawa couplings and masses to be fixed in terms of the light 
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Figure 4: Scatter plot of the muon EDM against the left-smuon mass for the ansatz Hi, 
taking the SU(2) gaugino mass to be 200 GeV, A = — 3m , /i > ; and tan/3 = 10. We 
assume the normal hierarchy for the light neutrino mass spectrum in (a) and the inverted 
hierarchy in (b). 



normal hierarchy inverted hierarchy 




m~[GeV] m~[GeV] 

Figure 5: Scatter plot of the electron EDM against the lefts electron mass for the ansatz H 2 . 
Other input parameters are the same as in Fig. 0. 
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neutrino parameters and a Hermitian parameter matrix H. In the minimal supersymmetric 
version of the seesaw model, the matrix H can be related directly to low-energy physical 
observables. As a result, our parametrization is particularly suitable for comprehensive 
studies of the charged LFV processes and EDMs in supersymmetric models. 

As applications, we have studied the LFV r decays and the EDMs of the muon and 
electron in the minimal supersymmetric seesaw model. It is found that Br(r — > n{e)^f) 
could exceed the present experimental bounds, even when Br(/i — > ej) is suppressed much 
below the current limit. This implies that B factories have a possibility of discovering LFV 
r decays, since they may reach sensitivities Br(r — > /z(e)7)~ 10 -8 . The LHC may have a 
similar sensitivity, and a super B factory may reach the level 10 -9 for the same processes [22|. 
The LHC may also have a good sensitivity to Br(r — > fi(e)£ + £~). We show that Br(r — > fiee) 
and Br(r — > 3e) are about five times larger than Br(r — > 3/i) and Br(r — > e2yu), due to the 
larger phase space, and can reach ~ 10~ 8 and ~ 10~ 9 , respectively, from the experimental 
bounds on Br(r — > /i(e)7). Finally, in our random samples the EDMs of muon and electron 
can attain lO^ 27-28 -^ cm and 10~^ 29 ~ 30 ^e cm, respectively, while their perturbative bounds 
in this model are ~ 10 _25 e cm and ~ 10 _27 e cm. The electron EDM, in particular, may be 



accessible to experiment p5 ]. 

We have restricted our discussion to the supersymmetric seesaw model in this paper. 
However, our framework is also suitable for studying supersymmetric GUT models with 
heavy singlet neutrinos. In these models, K and B physics are also interesting, because 
the right-handed squarks have flavour-violating masses, as a result of quark-lepton unifica- 
tion [p6fl . Since the right-handed sleptons have LFV masses, the relation between the LFV 
and our parametrization may be quite complicated However, if the prediction for the 
LFV masses for right-handed sleptons is used, our parametrization is applicable. 
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